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We provide a theoretical description of the polarization pattern and phase singularities experi- 
mentally evidenced recently in a condensate of indirect excitons [H. High et al., Nature 483, 584-588 
(2012)]. We show that the averaging of the electron and hole orbital motion leads to a comparable 
spin-orbit interaction for both type of carriers. We demonstrate that the interplay between a radial 
coherent flux of bright indirect excitons and the Dresselhaus spin-orbit interaction results in the for- 
mation of spin domains and of topological defects similar to Skyrmions. We reproduce qualitatively 
all the features of the experimental data and obtain polarization pattern as in the optical spin Hall 
effect despite the different symmetry of the spin-orbit interactions. 
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Introduction. The recent implementation of synthetic 
gauge fields acting on spinor atomic Bose Einstein con- 
densates [H has opened a new direction of research Q . 
Indeed, the coupling between the spin degree of freedom 
and the motion of bosons, linear in fc, completely mod- 
ifies the nature of the condensed states. The bosons do 
not accumulate any more in the single k = state, and 
the condensate can demonstrate spin-dependent modula- 
tions, such as stripes In confined systems, skyrmions 
and half vortices are expected to appear, and not as the 
excitations of the condensate, but directly in the ground 
state Q ■ The exact nature of these original equilibrium 
phases depends both on the type of spin-orbit coupling 
and on the inter-particle interactions. 

Indirect excitons in coupled quantum wells arc bosonic 
particles which continue to attract a strong interest since 
the Lozovik and Yudson's proposal to use a conden- 
sate of such particles towards a new type of supercon- 
ductivity [Hi- The spin structure of indirect excitons is 
especially interesting and unique since, differently from 
bosonic atoms, they possess an even number of allowed 
spin projections on the sample's growth axis. Two spin 
projections ±1 can couple to light to form bright states, 
while the ±2 components are bound to be dark. Further- 
more, both the electrons and the holes are affected by 
the Rashba and Dresselhaus spin-orbit interaction (SOI) 
which leads to a variety of accessible phases, both in the 
homogeneous Q and in the confined systems 0- 

However, the experimental phenomenology of indirect 
exciton condensation is complex and still not completely 
understood. Their principal advantage with respect to 
direct excitons or exciton-polaritons is the reduced over- 
lapping of their wave-functions that infers the quasi- 
particles a long lifetime in the range of nanoseconds. This 
feature is expected to optimize the process of thermal- 
ization of excitons which in turn should make the BEG 
formation easier. On the other hand, the critical conden- 
sation temperature is of about one Kelvin, a temperature 



at which an electron or an exciton gas is completely lo- 
calized by the structural disorder present in any semicon- 
ductor structure. Moreover, only a very tiny fraction of 
the exciton distribution can be probed by optical means. 
The dark states cannot be detected, and the bright ones 
are allowed to emit only within the light cone which rep- 
resents a very narrow part of the reciprocal space. De- 
spite these apparent difficulties, some interesting experi- 
mental results have been recently obtained in Ref. and 
can be summarized as follows. Far away from the lo- 
cal pumping area, strong maxima in the bright exciton 
distribution are observed. They were attributed to com- 
plex processes of carrier diffusion. These localised bright 
spots (LBS) were found to generate a radial flow of bright 
excitons characterized by: 1) a spatial coherence, with a 
coherence length of the order of 10 microns; 2) a specific 
linear polarization pattern over the same length scale; 
and 3) the presence of phase singularities. This specific 
linear polarization pattern has been attributed to the ra- 
dial cloud expansion from localized sources (See supple- 
mental material of Ref. 0] ) . The wave vector distribution 
of the electron and the hole wavefunctions forming the 
excitonic wave packet was neglected. This approach al- 
lowed to reproduce the polarization pattern, but not the 
presence of phase singularities observed. 

In this Letter, we consider a different and more com- 
plete approach, allowing to qualitatively reproduce all 
the experimental features. First, we take into account 
the internal structure of excitons and show how the rel- 
ative motion of electron and hole affects the SOI. Espe- 
cially, the hole SOI is found to become comparable to 
that of the electron. Second, we assume that a bright in- 
direct exciton condensate with a well defined spin state 
is locally formed at the center of the LBSjas assumed in 
Ref. [1] . Then, unlike the approach of Ref. Q , we consider 
the coherent expansion of the resonantly created exciton 
cloud by direct numerical solution of the Schroedinger 
equation, in the presence of the Dresselhaus SOI. The 
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renormalization of the dispersion induces a radial flow 
of excitons outwards from the pump spot. The repulsive 
exciton-exciton interactions can also contribute to this ef- 
fect, but we neglect them in the present work. Since the 
typical time scale of the scattering of indirect excitons on 
phonons is of about a few ns @ , the ballistic propagation 
length is expected to be of the order of 10-20 microns. 
Within this length scale, we reproduce both the polar- 
ization pattern, and the presence of phase singularities 
of the wave function components, which are associated 
with the formation of Skyrmions. These topological de- 
fects appear thanks to the interplay between the radial 
flow and the SOI, as recently shown theoretically for 
cavity polaritons flowing in a TE-TM effective magnetic 
field. Inspired by the experimental results, we propose a 
configuration leading to the onset of circular polarization 
domains, fully equivalent to the one observed in the op- 
tical spin- Hall effect This happens despite the fact 



that the wavevector dependence of the effective magnetic 
field of the Dresselhaus SOI is completely different from 
the one given by the TE-TM splitting in the microcavi- 
ties. 

The model. Indirect excitons are composite bosons 
formed from electrons and heavy holes separated in 
space. Condensates of indirect excitons were observed 
in systems of double quantum wells 12, l^j, where the 
electrons and the holes were localized in different layers, 
and in wide single quantum wells under applied or inter- 
nal electric fields [1J|. The projection of the electron's 
spin (angular momentum of heavy hole) on the growth 
axis can take two values ±1/2 (±3/2), so the total an- 
gular momentum of exciton can take four values ±2 and 
±1. The bright |±1) states are coupled with cr* -polarized 
light, while the states |±2) are dark. Their radiative re- 
combination is forbidden by selection rules. 

In planar structures these four states can be cou- 
pled via the inter-exciton interactions, the application of 
an external magnetic field, or by spin-orbit interactions 
(SOI) of carriers originating from the violation of the in- 
version symmetry. The interactions of the first type are 
density-dependent. There are two kinds of spin-orbit cou- 
pling for carriers: Dresselhaus and Rashba [3l SOIs. 
Both interactions act as effective in-planc fc-dependent 
magnetic fields coupling different carrier spin states |17|. 
However, in the sample used in the experiments of ref- 
erence [3, 0, , the Dresselhaus SOI has been found to 
be much larger than the Rashba SOI, and thus the latter 
shall be neglected in the following. 

For a gas of free 2D electrons, Dresselhaus SOI is linear 
versus the electron wave-vector and the corresponding 
term of the Hamiltonian in the basis (+1/2, —1/2)-'" takes 
the form: 



(1) 



where /3e is a Dresselhaus interaction constant. Here and 
after we use the following notations: kl = kj ± ik2, kZ „ 



are x- and y-components of wave- vector and 7 =c, h 
or X for electrons, holes or excitons respectively. a± = 
Gx ± icTy when ax.y are Pauli matrices. 

For free heavy holes in 2D, SOI terms are cubic versus 
the wave-vector fc, and in the basis (+3/2, —3/2)-^" have 
the following form [iq} : 



= — /3/l((T_|-fc_fc_|_fc_ +fT_fc+fc_fc+) 



(2a) 



(2b) 



Here au and (3h are interaction constants. Due to the 
cubic dependence on fc, SOI of holes was considered for 
a long time as a negligible mechanism of hole spin re- 
laxation [i^l- In latest studies 21 1, a splitting linear in 
k between heavy hole states has been predicted consid- 
ering the interactions of heavy-holes with remote energy 
states in the valence band. However, to our knowledge, 
nobody has yet considered how the Hamiltonian ([2]) is 
transformed for a hole bound inside an exciton. To do 
this, we decouple the relative motion of electron and holes 
from the motion of the excitonic center of mass: 



*(fc^,g) =«'c.„,(fc-^)*rcl(g). 



(3) 



Here q is large and rapidly varying wave- vector, describ- 
ing e-h relative motion and k^ relies to hole wave-vector 
as k^ = Vhk^ where v^^^h = 'm'e,h/{^e + ^h) is the ratio 
of the electron (hole) mass to the exciton mass. Then av- 
eraging Hamiltonian ([2]) over the relative motion ground 
state gives rise to Dresselhauss SOI linear in fc'* term for 
holes: 

H^ = I \^rcl{q)\'H^{k'' + q)d^q = h{<J+k1+a^k1), 

(4) 

while the Rashba term under similar averaging remains 
cubic in k^ and can be anyway safely neglected. Here 
I3h — —2/3/1 {q'^) is the effective Dresselhaus constant for 
holes. Note that (g^) ~ a^^ » (k'^)'^, where as is the 
excitonic Bohr radius. Analogous procedure for electrons 
does not change the form of ([T]) but renormalizes electron 
interaction constants. Finally, one may write a SOI term 
for both electron and hole bound into exciton as: 



= /3e,hi^eAcr+k^ + a-kf). 



(5) 



Here /3e,h are the effective Dresselhaus interaction con- 
stants of electrons and holes. 

In addition to the SOIs, fc-independent energy splitting 
between linearly polarized states, e.g. parallel and per- 
pendicular to crystallographic axis can occur, and can 
be described as /c-independent constant effective mag- 
netic fields. Finally, we can write the Hamiltonian for 
indirect excitons on the basis of the four spin states 
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(6) 
are the 



parabolic dispersions of the bare indirect exciton states, 
6b,d give the energy splittings between linearly polar- 
ized bright (dark) states. We take the exciton energy 
at k-^ = as the zero point. Putting 5b = 5d = 
the Hamiltonian diagonalization is straightforward and 
yields the following set of isotropic eigenmodes: 







/ Eo 






+ VePe) k^ 


\ 


En (fc-^) 




Eq 




+ {vhPh 


- yePe) fc^ 




Eiii (k^) 




Eq 




- {i^hPh 






V Eiv (fc^) ) 




\Eo 




- {vhPh 


+ VePe) k^ 


J 



(7) 



and the corresponding cigcnstates: 
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(8) 



Here </> is the polar angle in reciprocal space. The combi- 
nations of EqsIH correspond to "linear polarizations" for 
the bright (rows 2 and 3) part given the usual identities 
tjjj. = tpx i iipY- While dark components are isotropic, 
the bright components have a linear polarization that is 
0-dependent and the polarization changes from X to y 
when (f) is changed by 7r/2. This peculiarity is analogue to 
the exciton-polariton case in the presence of the so-called 
TE-TM splitting that gives birth to the optical spin-Hall 



effect I 111 . |22| . |23| . One can therefore expect the for- 
mation of the same polarization pattern from the bright 
exciton states under proper excitation. 

The isotropic dispersion branches given by the Eq.© 
are plotted on the Fig H] in the case 6b = Sd = 0- One can 
clearly see the energy splittings between the branches 
coming from the Dresselhaus SOI contribution. Ex- 
perimentally, one can detect only the emission coming 
from bright excitons with small wave-vectors inside the 
light cone (2.6 x 10^m~^). In this region the dispersion 
branches are linear. Interestingly, the calculated ground 
state of excitons has a significant wave vector completely 
out of the light cone. Thus a condensate in the ground 
state, even possessing some "bright" exciton component, 
is expected to be completely dark. It is moreover de- 
generate, and may demonstrate rich phenomenology of 
topological defects [24| . The Liouville equation allows us 




FIG. 1: (Color online) Isotropic dispersion branches of the 
indirect exciton eigen modes [see Eqs.((7])]. 



to derive the reciprocal space dynamics of the spinor field 

+ 1: 



* = (V'+2,'0+i,'0-i,'0-2)^ with Eq.® yielding 
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where we have phenomenologically introduced the exci- 
ton decay with a lifetime r and a local source P(k-^) = 
Ap (k-^) (5 (cj — cjp) with frequency up (in our calcula- 
tions we took hujp ~ 2^eV ) acting on each compo- 
nent, where Ap(k-'^) are 2D Gaussians. The forma- 
tion of the condensate of indirect excitons coming from 
a hot reservoir and the process of their relaxation to- 
wards ground state are complex and need further inves- 
tigations. Within this simplified model we consider a 
narrow pump spot (broad in reciprocal space) exciting 
laterally the dispersion branches and assume a ballistic 
propagation of the exciton. As said above, the ballistic 
propagation time is of order of a few ns, during which the 
evolution of the created cloud is expected to be coherent. 

In the idealized case of a Dirac delta source in real 
space, the excitation of an eigenstate can be found an- 
alytically. For example, the dynamical equation for 
ipi{k^) reads: 



ih 



^^PI{k^,t) _ h^ 



ik^'fi^j (k^,t) 



dt 2m' 
+iA(3k^iji (k^,t) + Ape"^''* 



and the stationary radial solutions are found writing 
^pi{k^ ,t) = tpi {k^ ) exp{iujpt) which yields: 
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the Green's function of the problem reads 
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whose Fourier transform gives 
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which is nothing but Tjji{r) for the delta source. K1.2 are 
the poles of Ea.([TT|) and Eas. (|12ll3p require that ki > 
and K2 < which means that hup > 0. The solutions for 
the other eigenstates are found with a similar procedure 
and one can therefore construct any combination using 
the eigenstates ([8]). 

For a numerical simulation, we first consider a linearly 
polarised {P±2 ~ 0, P±i ^ 0) bright cxciton state as ini- 
tial condition. Here and after we took equal Dresselhaus 
constants for electrons and holes (3e,h — QfJ.eV ■ fim, and 
the masses rUe ~ 0.07mo, m/i = 0.5?tio, where mo is the 
free electron mass. The figure [2j a, b,c) shows the station- 
ary real space images of a linear, circular polarization 
degree and phase of bright states. One can observe the 
formation of a Skyrmion lattice associated with the for- 
mation of the spin domains. The situation is analogous 
with the polaritonic spin Hall effect recently analyzed 
theoretically [l3l- This similarity is expectable, since the 
effective field acting on the bright states is exactly equiv- 
alent to the case of exciton polaritons in the presence of 
the so-called TE-TM splitting. FigH (c) shows the den- 
sity profile of bright states along the diagonal line [white 
dashed line on Fig[2fb)]. Densities of a"'" and a~ states 
are oscillating while the total density decays in space with 
r. Phase structure of tr^-polarized bright state is plot- 
ted on Fig[2I^d). However, the linear polarization pattern 
that we observe (Spolarization domains), differs from the 
one measured in [8[ (4 domains). 

In order to reproduce the experiment, we therefore 
consider a different initial spin state for the conden- 
sate (Fig. 3). We consider first a condensate of dark 
states with a slight asymmetry between dark components 
{P_2/P+2 = 0.9, P+i = P-i = 0). Then, in order to mix 
the circularly polarized bright states, we introduce an 
additional constant splitting between the linearly polar- 
ized bright states along and perpendicular to the main 
crystallographic axis of the sample 5h = l^eV . 

Fig. [3] (a,b) shows the spatial distribution of linear 
and circular polarization degree, while FiglSJd) shows 
the phase of the X-component. Both circular and lin- 
ear polarization structures present 4 domains and are 
stretched along the x axis (because of 5b) as it was ob- 
served in experiment 0. Interestingly, we also observe 
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FIG. 2: (Color online) (a) Degree of linear polarization of 
bright states, (b) Degree of circular polarization of bright 
states, (c) Density profile for a'^- {nb+, red dashed line),cr^- 
{ub- , blue dashed line) polarized bright states and total den- 
sity (n^"*, black solid line) of bright states along the diagonal 
line y = —x (white dashed line in (b)). (d) Phase of a^- 
polarized bright component. 
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FIG. 3: (Color online) (a) Degree of linear polarization of 
bright states. (b)Degree of circular polarization of bright 
states, (c) Total density of bright states, (d) Phase of x- 
component of bright states. Red crosses point the phase sin- 
gularities in s-component 



pairs of phase singularities (red crosses), situated sym- 
metrically with respect to the exciton source. This phase 
singularity is accompanied by a density dip only in the 
Y component. This topological defect is therefore sim- 
ilar to a Skyrmion but in the linear polarization basis. 
Our approach shows that the appearance of phase singu- 
larities is a general feature of radial flows of particles in 
the presence of coupling between the spin and motional 
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degree of freedom. 

Conclusions. We have shown that an expanding cloud 
of indirect excitons can form polarization textures and 
phase singularities thanks to the interplay between the 
flow and spin-orbit interactions. Their structure strongly 
depends on the polarization of initial state, so chang- 
ing the excitation conditions one can strongly modify the 
number of polarization domains and qualitatively repro- 
duce experimental data. Additionally, we were able to 
mimic the optical spin-Hall effect for bright excitons. 
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